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A macroscopic  model  of  heterogeneous  nucleatlon  is  used  for  a theo- 
retical study  of  condensation  of  water- vapour/carrier-gas  mixtures  in  a non- 
equUihrlum  nonstationary  rarefaction  wave  generated  in  a shock  tube.  The 
results  are  conpared  with  those  from  homogeneous  nucleatlon.  Nucleatlon  is 
assumed  to  take  place  heterogeneovisly  on  idealized  smooth,  spherical  solid 
particles  of  Aitken  nuclei,  which  are  chemically  and  electrically  inert.  In 
the  processes  of  heterogeneous  condensation,  the  controlling  factors  are  the  i 

size-distribution  of  nuclei,  the  concentration  of  monomers  on  the  surface  of  i 

the  substrate  and  the  contact  angle  of  embryos.  Of  these  factors  the  most 
dctninant  is  the  contact  angle,  which  can  reduce  greatly  the  activation  energy  j 

of  nucleatlon.  Heterogeneous  condensation  results  in  less  siq)ercooling  of  the 
mixture  bixt  a f sister  approach  to  the  eqiilllbrlm  state.  By  choosing  a suitable 
value  for  the  contact  angle,  the  nianerical  results  can  be  made  to  fit  the 
experimental  data.  Although  this  is  not  entirely  satisfactory,  it  is  probably 
preferable  to  changing  the  value  of  surface  tension  as  in  the  homogeneous 
nucleatlon  case  in  order  to  obtain  agreement  with  experimental  results. 
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space  coordinate 
y Lagrangian  coordinate 

O condensation  coefficient 

nondimensional  adsorption  free  energy,  Eq.  2h 
7 specific  heat  ratio 

0 enibryo-to-substrate  contact  angle 

X nondimensional  latent  heat  of  evaporation 

p Lax  dissipation  constant 

V adatom  vibrational  frequency 
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p density 
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1.  IMTRQDUCTION 


Problems  of  two-phase  and  tvo-conponent  systems  have  great  practical 
Inportance  in  various  engineering  fields.  Resecurches  on  the  nucleation  and 
growth  by  iconcjnsation  and  evj^joration  of  droplets  have  been  conducted  over  a 
long  period  of  time  and  many  aspects  of  these  processes  can  be  regarded  as 
understood.  Nevertheless,  there  still  remains  a number  of  unsolved  problems 
such  as  phase  transitions  under  the  existence  of  foreign  particles. 

Homogeneous  nucleation  and  condensation  have  been  stvidled  theoretically 
and  experimentally  by  a nuidser  of  investigators  within  the  framework  of  one-  or 
two-dimensional  siper sonic-nozzle  flows  and  nonstationary  rarefaction  flows  in 
a shock  tube.  Recently,  Sislian  (Ref.  l)  made  a detailed  numerical  study  of 
homogeneous  nucleation  and  condensation  of  water  vapour  with  or  without  a 
carrier  gas  in  the  nonstaticnary  rarefaction  wave  generated  in  a shock  tube  and 
predicted  the  effects  of  condensation  on  the  flow  variables. 

Unless  a siper saturated  vapour  is  specially  and  expensively  treated, 
it  is  likely  to  contain  a large  nusber  of  small  particles  (Altken  nuclei), 
which  can  act  as  nuclei  of  spontaneous  condensation.  Depending  upon  the 
cooling  rate  and  the  number  of  nuclei,  the  resulting  phase  change  of  the  svper- 
satvurated  vapo\ir  can  take  place  at  considerably  lower  svper  saturation  theui 
predicted  by  homogeneous  nucleation  alone. 

Since  nucleation  on  foreign  particles  Involves  additional  degrees  of 
freedom,  the  anelysls  of  heterogeneous  nucleation  is  much  more  conplicated  than 
the  homogeneous  case.  The  nucleation  rate  will  depend  on  the  natvire  of  the 
nucleus  surface,  including  its  geometry  and  physical  and  chemical  properties. 

It  also  depends  on  the  condition  of  the  surfeice  and  the  state  of  the  vapour. 
Nevertheless,  there  are  sany  similarities  between  homogeneous  and  heterogeneous 
nucleation.  Due  to  these  featvires,  heterogeneous  nucleation  has  been  studied 
from  both  macroscopic  and  microscopic  points  of  view,  such  as  the  capillarity 
model  (Refs.  2,  6)  and  the  small-cluster  model  (Ref.  3)*  The  latter  requires 
fairly  microscopic  data  of  physical  and  chemical  properties  of  the  nucleus 
surface  which  are  not  readily  available. 

In  the  present  report,  by  using  a macroscopic  model  of  heterogeneous 
nucleation,  numerical  studies  are  made  of  the  condensation  of  water  vapo\ir  in 
the  nonstationary  rarefaction  wave  generated  in  a shock  tube  and  the  results 
are  conpared  with  those  obtedned  from  homogeneous  nucleation  (Ref.  l).  Nuclea- 
tion is  assumed  to  take  place  heterogeneously  on  idealized  smooth,  spherical 
particles  of  Altken  nuclei,  which  are  chemically  and  electrically  inert.  The 
flow  variables  in  a shock  tube  are  assumed  to  have  the  same  initial  conditions 
used  in  Ref.  1.  The  new  flow  properties  are  then  c€d ciliated  using  the  same  Lax 
method  enployed  in  Ref.  1. 

2.  HETEROGENEOUS  NUCLEATION  AND  CONDENSATION 

2.1  Condensation  Nuclei 

As  a source  of  commercial  nitrogen  gas,  atmospheric  air  is  widely  used. 
The  air  has  sufficient  quantities  of  nuclei  produced  by  physical,  chemical  and 
mechanical  processes  in  nature.  According  to  cloud  physics  terminology  (Ref.  4), 
these  particles  eu'e  classified  by  dimension  as,  Altken  nuclei  for  particles  with 
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radii  below  O.ln;  large  nuclei  for  those  with  radii  0.1  ~ In;  and  gigantic  or 
giant  nuclei  for  those  with  radii  in  excess  of  I41.  Large  and  gigantic  nuclei 
usually  conteiin  a large  quan\.::uy  of  soluble  substances,  primarily  ansnonium 
sulfate  and  sodium  chloride.  They  are  the  most  active  of  nuclei  in  cloud 
formation.  The  majority  of  Aitken  nuclei  are  too  !=™a.n  to  become  active  in 
meteorological  condensation  provided  they  €ire  not  electrically  charged. 

Even  cleaned  air  with  careful  filtration  might  carry  a considerable 
nunber  of  these  condensation  nuclei.  In  industriaO.  processes,  the  gases  are 
passed  through  various  pieces  of  equipment,  heat  exchangers,  condensers,  liq\iid- 
sprayed  colianns,  contact  apparatus  with  cateQ.yst  layers  and  filters.  The  number 
density  of  pairticles  and  their  physical  and  chemical  properties  thus  depend  not 
only  on  those  in  the  initial  gas  but  also  on  the  processing  of  the  gas.  By  taking 
special  care  of  physical  filtration  and  electrical  neutralization,  large  and 
gigantic  particles  and  electrically  chaurged  particles  can  be  excluded  from  the 
gas.  Consequently,  the  processed  gas  can  be  assumed  to  contain  only  Aitken  nuclei 
that  are  chemically  and  electrically  inert. 

2.1.1  Number  Density  and  Size  Distribution 

Concentration  and  size  distribution  of  atmospheric  nuclei  vauy  greatly 
with  locality.  The  particles  in  greatest  concentration  are  the  Aitken  nuclei, 
which  are  counted  in  the  conventional  Aitken  nucleus  coiinter,  ranging  in  nvunber 
density  of  10^  to  10^  (pairticles  per  cm^).  A representation  of  the  average  size- 
distribution_in  typical  atmospheric  air  is  shown  in  Fig.  1 (kef.  5).  Two  peaks 
below  3 X 10  ® cm  are  due  to  more  aggregations  of  air  and  water  molecules  on 
ions,  and  are  not  foreign  particles.  Hence,  the  gas  treated  with  physical  fil- 
tration and  electrical  neutralization  can  be  assximed  to  have  such  a size  distri- 
bution of  nuclei  given  by  curves  A,  B or  C. 

In  order  to  describe  the  number  density  and  size  distribution,  a 
distribution  function  with  respect  to  the  droplet  radius,  R,  sho\ild  be  defined. 

The  distribution  function  N(R)  multiplied  by  dR  characterizes  the  number  of 
particles  lying  in  the  size  range  from  R to  R + dR, 


dN  = N(R)dR 


The  overall  particle  concentration  is  then  given  by 


N > = J'  N(R)dR 


The  ass\aned  distribution  functions  in  Fig.  1 can  be  expressed  in  a form  of 


dN 

d?i^ 


where,  Nq  and  Rq  are  given  in  Table  1,  with  other  size  average  properties. 
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2.1.2  Shape  and  Size 


Particles  with  a radius  less  than  lO”  cm  (lOO  S)  are  highly  dispersed 
and  in  vigorous  Brownian  motion  so  that  the  shape  of  the  particles  is  not  at 
all  spherical  due  to  coagulation  and  sintering  phenomena.  It  may,  however,  be 
possible  to  classify  them  by  some  effective  particle  size  such  as  an  average 
diameter  or  radius,  although  small  pits,  steps,  cracks  and  cavities  in  the  surface 
are  the  most  ing)ortant  in5)erf actions  from  a physical  point  of  view  of  heterogeneous 
nucleation.  Nevertheless,  for  sinplicity,  all  particles  are  assvnned  to  have  a 
smooth,  spherical  shape  of  effective  radius  R,  with  a distribution  function  given 
by  Eq.  3,  without  any  further  coagulation  or  separation. 

2.1.3  Surface  Physics  and  Chemistry 

The  main  sources  of  Aitken  nuclei  are  smoke  and  vapour  from  fires  and 
industries,  dust  from  the  lands,  salts  from  the  oceans  and  particulate  products 
from  chemical  reactions.  The  chemical  and  physical  properties  of  these  nuclei 
are  extremely  varied.  This  arises  because  the  fonnation  of  the  nuclei  is  due 
not  only  to  a number  of  their  originating  processes  but  also  to  subsequent 
coalescence  of  particles,  adsoiption  of  gases,  and  chemical  reactions.  Hence, 
one  may  actually  have  reference  only  to  average  properties  of  particles. 

If  the  nucleus  consists  of  a mixture  of  soluble  and  insoluble  particles, 
then  nucleation  is  dominated  by  the  soluble  particles  due  to  their  smaller  acti- 
vation energy.  The  gases  treated  with  physical  and  chemical  processes  can, 
however,  contain  a negligible  quantity' of  soluble  particles.  Here,  all  condensa- 
tion nuclei  are  assumed  to  be  chemically  inert  and  insoluble  in  water  vapour. 

2.2  Heterogeneous  Nucleation 

From  the  macroscopic  point  of  view,  heterogeneous  nucleation  follows 
steps  very  similar  to  homogeneous  nucleation,  although  the  fomer  occurs  on  a 
surface  immersed  in  a supersaturated  vapour.  Bsvereil  models  have  been  proposed 
for  heterogeneous  nucleation  (Refs.  6,  7).  'According  to  these  models,  nucleation 
takes  place  as  a result  of  growth  beyond  a critical  size  of  embryos  from  a thermo- 
dynamic eq\ailibri\mi  distribution  colliding  with  monomers  or  higher-order  clusters. 
In  addition,  embryos  can  also  form  from ' adsorbed  monomers  or  higher-order  clusters 
in  one  or  more  layers  on  the  solid  surface. 

The  nucleation  rate  per  unit  siorface  area  of  particle  per  unit  time  can 
be  basically  expressed  as, 


n = Z • w • n(i*)  (k) 

where,  n(i*)  is  the  concentration  of  embryo  of  size  i *,  on  the  substrate;  w, 
the  in5)ingement  rate  of  monomers  on  the  eribryoj  and  Z,  the  nonequilibrium  correc- 
tion factor. 

2.2.1  Nucleation  Mechanisms 


Consider  monomers  or  higher-order  clusters  of  a sx^qjersaturated  vapour 
in5)inglng  on  a surface  of  a solid  particle  (R),  to  make  a cap-shaped  embryo  (r), 
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as  illustrated  in  Fig.  2.  For  heterogeneous  nucleation,  there  are  two  types 
of  monomer  ispingements  on  embryos;  the  direct  ispingement  of  monomers 
from  the  vapour  and  the  inpingement  by  desorption  of  adsorbed  monomers. 

(i)  Direct  inpingement  of  monomers  from  vapoiar  (Ref,  7) 

According  to  kinetic  theory,  for  a Maxwellian  velocity  distribution, 
the  monomer  flux  J(l)  is  given  by. 


where,  pv  is  the  vapovir  pressure,  T the  tenperature,  m the  mass  of  one  molecule 
and  k the  Boltzmann  constant.  The  inpingement  rate  of  monomers  on  the  entoryo 
is  then, 

V = J(l)a*  (6) 

where,  ai^  is  the  surface  area  of  the  enibryo. 

The  concentration  of  embryos  of  the  size  i*  (number  of  monomeisi*)  can 
be  epproximately  given  by, 

n(i*)  = ng(l)  exp  j-  (7) 


where,  ng(l)  is  the  concentration  of  adsorbed  monomers  on  the  substrate  and 
2!!G*(i)  is  the  Gibbs  free  energy  of  formation  of  an  embryo  of  size  i*,  from  the 
vapour  phase  monomer  of  concentration  ng(l) . 

The  concentration  of  monomers  adsorbed  on  the  substrate  can  usxially 
be  written  in  terms  of  the  inpinging  vapour  flux  J(l)  and  the  mean  manomer  stay 
time  on  the  substrate  t as, 

n (l)  = a J(1)t  (8) 

s a 

where,  a is  the  adsorption  coefficient  and  t is  given  by, 

T =1  V exp  (^-  (9) 

where,  AG^  is  the  desorption  free  energy  of  monomers  and  v is  the  vibrational 
frequency  of  the  adsorbed  atom. 

The  nucleation  rate  per  unit  surface  area  of  a nucleus  particle  by 
direct  inpingement  of  monomers  from  the  vapour  is  then. 
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(ii)  Inpingement  of  adsorbed  monomers  on  the  substrate  (Refs.  7,  8) 

The  rate  at  'which  adsorbed  monomers  Is^inge  upon  embryos  on  the  sub- 
strate is  given  by  (Ref.  8) , 

w = J(1)t  V exp  (11) 

where,  ££  is  the  embryo  periphery,  is  mean  siurface  diffusion  distance  and 
AGg  the  free  energy  of  activation  for  surface  diffusion. 

By  sTibstituting  Eqs.  5,  7,  8,  9 and-  11  into  Eq.  4,  the  nucleation  rate 
can  be  •written  as. 
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Using  the  reasonable  assumption  that  Ajg  « and  ii  ig  *»  both 
nucleation  rate  eqviations  can  be  expressed  in  the  same  form  as. 


- nAG^ 


where,  n = 1 for  the  direct  impingement  of  monomers  from  vapour  and  n = 2 for 
the  impingement  of  desorbed  monomers.  Thvis,  the  difference  between  nucleation 
mechanisms  of  direct  and  indirect  impingements  of  monomers  can  be  included  in 
the  desorption  free  energy  from  a macroscopic  point  of  view  (nAO^  -» AGd) . In 
this  sense,  the  nucleation  rate  per  unit  area  of  substrate  csui  be  expressed  as. 


AG*  - AG, 


(iii)  Noneqviilibrium  correction  factor 

Because  the  critical  nuclei  are  steadily  being  cansvuned  to  produce 
droplets,  the  actual  concentration  of  emibryos  and  critical  nuclei  are  smaller 
than  the  equilibrium  values.  The  factor  which  corrects  this  embryo  depletion 
is  given  by. 


for  a cap- shaped  embryo  (Ref.  7)>  although  it  is  not  ordineurily  a significant 
quantity  in  the  nucleation  rate. 


(iv)  The  critical  free  energy 
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For  embryos  of  small  size,  the  partition  functions  for  monomers  and 
clusters  including  the  translational,  vibrational,  rotational  and  configura- 
tional contributions  should  be  ev6^.uated  to  describe  the  free  energy  of  embryo 
formation.  However,  because  of  its  sinplicity  and  mathematical  tractability, 
a macroscopic  approach  is  often  vised.  By  this  method,  the  free  energy  of 
embryo  formation  can  be  expressed  as  the  svnn  of  contributions  from  the  forma- 
tion of  the  solid- liqviid  interface,  the  liquid-vapovir  interface,  and  the  bvilk: 


where,  is  the  free-energy  difference  per  unit  volvime  of  liqviid  phase  from 
v^our  to  liqviid  phase,  the  surface  free  energy  of  the  interface  between 
phases  i and  d,  v^,  the  volume  of  liquid,  a,.,  the  svirface  area  of  liquid-vapovir 
interface,  ag,  the  surface  area  of  liquid- substrate  interface  and  subscripts 
V,  i and  s refer  to  vapour,  liqviid  and  substrate,  respectively.  From  geometri- 
cal consideration,  v£,  a^-  and  ag  are  written  as, 

v^  = ^ nr^  s.  cose  ^ 

= Ftln-2  cose  ^ (17) 

where,  a,  by  and  bg  are  a function  of  r/R  and  the  contact  angle  0,  given  in 
Appendix  A,  and  in  Fig.  3. 

By  using  the  usual  definition  of  the  contact  angle 


a . - a 

cose  = SZ 


a . - a 
s£  sv 


(■"  ' ■'ilv) 


(18) 


and  substituting  Eqs.  17  and  l8  into  Eq.  l6,  the  free  energy  of  fonnation  of  an 
embi^o  of  radius  r , on  a nucleus  of  radius  R,  can  be  expressed  as  a function  of 
r,  R and  cos6.  Differentiating  /iG  with  respect  to  r to  evaluate  the  critical 
free  energy  AG*  gives  the  critical  radius  r*. 


at  which. 


6 


I 


The  critical  radivis  is  identical  with  that  for  homogeneous  nucleation.  The 
free-energy  difference  between  the  vapour  and  liquid  phases 


= - n^  kT  in 


kT  in  s 


where,  pg  is  the  saturation  pressure  of  the  vapour  phase,  n^,  the  nuntoer  of 
molecules  per  unit  volume  of  liquid  and  s,  the  supersaturation, 


SiJbstituting  Eqs.  19  and  20  into  Eq.  l6  gives  then  the  free  energy  of  formation 
of  a critical  embryo  (Ref.  9)> 

IStju®  ^ f t*  ^ 


IStju®  ^ f t* 

t£t*  = ' ' ■'  f { COS0  I 

3(n^  kT  ins.)^  > ^ ' 


where. 


f (|,  co.e)  - |[i  .{  c,  (f  . cose  )]■  + (£)  {s  . 3cg(l  Tj  cose) 
* c|  ( X - 1 cose  J}  + 3 (ijcose  I Cg  ( 1 - I cose  ) - i (23) 

The  factor  f(r/R,  cos6)  is  shown  in  Fig.  4.  It  means  that  lower  values  of 
contact  axigle  reduce  greatly  the  activation  energy  of  the  formation  qf  critical 
embryo  by  virtue  of  the  difference  of  surface  free  energy  between  vapovir-solid 
and  liquid-solid  Interfaces. 

Finally,  by  combining  Eqs.  l4,  15  and  20,  the  heterogeneoiu  nucleation 
rate  per  unit  area  of  nucleus  substrate  can  be  expressed  as, 

* ' {• 
where,  n^  mews  the  homogeneous  nucleation  rate  per  unit  volume  of  viq>our. 


n = 4T7r*^  " 

° sIzmlkT 


^v  f 47T 
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^ vV  27nn  ® ~ W 

o 


2.2,2  Nucleatlon  Rate 


A nucleus  particle  has  the  surface  area  Aq(*  kifR^)  at  the  initial  stage 
of  nucleation  and  the  area  cooes  to  be  wetted  by  condensed  droplets,  say  Ao  at 
time  t.  The  surface  area  of  a nucleus  on  'vdiich  embryos  can  nucleate  is. 


A^  = A - A 
e o s 


Since  the  size  distribution  function  of  nuclei  is  N(R)  and  the  nucleation  rate 
per  \init  surface  area  of  a nuclei  is  A,  the  total  nucleation  rate  per  unit  volume 
of  vapour  is  then. 


i =J^ A Ag  N(R)dR  = J^n(A^  - A^)  N(R)dR 


The  wetted  surface  area  A is  given  by, 

s 

Ag  = (^o  ■ -^s^T  (27) 

I ° 

where,  ag(t,r)  is  the  liquid-solid  interfacial  area  of  a condensed  droplet  idiich 
is  formed  at  time  t (<  t).  The  solution  of  the  integral  equation,  Eq.  27,  gives 
the  ^wetted  surface  area  of  a nucleus  Ag,  hence  the  total  nucleation  rate  A(A  - A ) 
or  I,  The  integral  equation  can  be  more  easily  solved  by  replacing  it  by  ° ® 

successive  differential  equations  (Appendix  B),  as  follows; 


= A(A  - A ) b + r B 

ot  ' o s'  • s os 


;nr-  = A(A  - a ) 2Trr*  c . + f C 
ot  ' o s'  . sb  o s 


" A ) 27T  c + r C 
ot  o s'  sc  o SI 


= A(A  - A ) 27T  c + r C 
ot  os  SCI  o S2 

• • • 

where,  Cg-jj,  Cg^,  Og^j^,,..  are  a function  of  r/R  and  cos0  given  by  Eq.  B-U  and 
ro  is  the  growth  rate  of  homogeneously  condensed  nuclei  given  by  Eq.  31  in  the 
following. 
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2.3  Condenaatlon  and  Droplet  Growth 


2,3,1  Droplet  Growth 

By  uelng  the  Hertz-Khudsen  equation  for  the  net  mass  flow  in  unit  time 
of  vapour  molecules  condensing  onto  a droplet  as  in  the  homogeneous  case  (Ref.  l) 
the  drqplet  growth  rate  can  he  approximately  e^Q>ressed  as. 


3F  (p/  = ‘v 


where,  pj  v/  is  the  laass  of  a condensed  droplet,  Sy>  its  liquid- vapour  surface 
area  and  a,  the  so-called  condensation  coefficient  eqvial  to  the  ratio  of  mole- 
cules sticklog  to  the  droplet  to  those  inplnglng  on  it.  The  above  relation 
with  Pjj  v|  « pjg  4TO®a/3  yields 


where,  fo  means  the  growth  rate  for  homogeneous  cases. 


° J2S5 


fC  1 + 


2 ,3 ,2  Condensate  Mass 


The  mass  of  embryos  produced  diuring  the  tine  r and  t + dr  is. 


J* ft(A^  - Ag)p^  v|  N(R)dRdT 


which,  grows  at  time  t vp  to. 


/” ■^/(^»'>’)N(R)dRdT, 


Thus,  the  overall  condensate  mass  per  unit  mass  of  vapovir  mixture  is  given  by. 


g = J”gp  N(R)dR 
o 

«R-  f r (*o  - 


and  the  rate  of  production  of  the  liquid  phase  by, 


^ = f 

dt  J 

o 


N(R)dR 


^ { I ''I  \ {I  - V ^ v}/'' 

o 

Further,  idien  homogeneous  nucleation  is  tahen  into  account  simultanecusly,  the 
rate  of  liquid-phase  production  can  be  written  as, 

i4 


■f[l 


(Ao  - 5 a* 


I K - V It  (5  ‘ fi) 


The  first  and  second  'Cerms  of  the  ri^t-hand  side  are  the  contribution  from 
homogeneous  nucleation. 

When  the  heterogeneous  nuclei  are  in  the  size  range  of  one  or  two 
orders  larger  theui  that  of  the  embryos  (R  ~ 10  cm)  and  the  values  of  contact 
angle  are  not  too  large-  (9  < 120°),  heterogeneous  nucleation  is  always  favoured 
compared  with  homogeneous  nucleation  because  of  its  smaller  activation  energy. 
For  the  same  sattiration,  the  heterogeneous  nucleation  rate,  n,  is  of  one  or  two 
orders  larger  than  the  homogeneous,  Hq  (see  Fig.  10).  Consequently,  in  the 
present  analysis,  only  the  heterogeneous  condensation  terms  are  considered  in 
the  rate  equation,  Eq.  34. 

The  rate  equation  can  be  numerically  solved  by  transforming  them  to 
successive  differential  equations  (Appendix  B),  as  follows: 


3T  = Pi 


)X  • 

-y  a + r 


i = I (A„  - A.)-*  % * ^o<= 
H - f °=  * "o<=l 


^ “ (A  - A ) c„  + r C_ 

dt  p ' o s'  Cl  o 2 
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where,  c^,  are  a function  of  r/R  and  cosS  given  by  Bq,.  B-7. 

Since  coefficients  a,  by.,  c's  in  Eq.  36  and  bg,  Cg's  in  Eq.  28  are 
a function  of  the  radius  of  nuclei  R,  the  total  mass  fraction  of  liquid  phase 
g should  be  obtedned  by  integrating  gp  with  respect  to  R,  which  is  extremely 
time-consuming  for  numerical  con^utstlon.  Hence,  it  may  be  sufficient  to 
consider  an  appropriate  approximation  for  the  integration  of  these  rate  equations, 
before  applying  them  to  the  problem  of  condensation  of  water  vapoiur  in  a shock 
tube. 

For  exaaple,  consider  the  system  of  water-vapour  and  air  mixtures, 
tdiich  contain  nuclei  particles  with  size  distributions  shown  in  Fig.  1 and 
have  the  initial  conditions  shown  in  Table  2.  The  pressure  of  the  system  is 
assumed  to  change  as 


^ = 0.425  + 0.525  e:q)(-  SOt^)  (37) 

which  is  nearly  similar  to  the  curve  at  -5  cm  from  the  diaphragm  in  Fig.  12a. 
By  combining  the  equations  of  energy  and  state,  Eqs.  53  and  4l,  the  rate 
equations,  Eqs.  28  and  36,  can  be  Integrated  numerically. 

In  Fig.  6 is  shown  the  variation  of  condensate  mass  fraction  with 
respect  to  time  corresponding  to  the  size  distributions  of  nuclei.  A,  B and  C 
in  Fig.  1.  In  order  to  examine  the  effect  of  size  distribution  of  nuclei,  the 
following  physical  properties  were  xued  in  Eq.  24: 


Rjj  < R®  > = 1,  = 34 


which  means  v = (l  ~ 10)»10^®/sec  smd  = 20  kcal/mol.  Eqviations  28  and  36 
were  integrated  vp  to  Cg  and  C,  respectively.  It  is  seen  that  because  of  the 
smaller  values  of  embryo  size  cospeured  i#lth  the  effective  radius  of  nuclei 
r*  « R,  the  size  distribution  has  little  Influence  on  the  result.  By  virtue 
of  this  fact,  coefficients  in  Eqs.  28  and  36,  a,  bv,  bs,  c's  and  Cg's  can  be 
sufficiently  approximated  with  their  values  in  the  limit  of  r/R  -»  0,  which 
are  given  by  Eqs.  A-4  and  B-12.  The  rate  equations  are  then  reduced  to  the 
following  ordinary  differential  equations: 


<^R  , 


H = - (1  - A ) r*  c + r C 
dt  p ' s'  b*»  o 


(39) 


dt 


c 

C«o 


U 


r 


dA 

IT  = “(1  - h'>  V * 'o  K 


cLB 

dT  = "(1  - Ag)  27rr* 

d5 

~ = ii(l  - A ) 27T  c _ 
dt  ' s'  SCOO 


(^) 


where 


®bo  = (1  - cose)^  ( I ■*■  ? ) 

V = I ^s«  = ¥ 


= 


W’ 


beo  2 + COS0 


°sboo=  (1  - 

K„  =1  I (1  + cose) 

and  f(r/R,  e)  involved  in  n is, 


^COO  ~ ^00 


^SCOO  ~ ^00  ^sb«o 


I + f cose 


)r 


= (1  - cose)^  ( I t °°se  ^ = 


Figures  7 to  10  show  the  effects  of  the  vibrational  frequency  adatoms 
(Rh),  the  desorption  free  energy  of  monomers  (^d)J  and  the  contact  angle  (e), 
upon  the  variations  of  ten5>erature,  supersaturation,  condensate  mass  fraction 
and  nucleation  rate  of  embryos  associated  with  the  change  in  pressure  given  by 
Eq.  37.  For  smaller  values  of  contact  angle,  heterogeneous  nucleation  occurs 
faster  at  lower  svq>ersat\iration  (Figs.  8 and  lO),  so  that  the  decrease  in 
temperature  is  smaller  due  to  rapid  heat  release  of  condensation  (Fig.  7) . 
Larger  values  of  the  factor  Rjj  < R^  > mean  larger  values  of  the  concentration 
of  adsorbed  monomers  on  the  nuclei,  or  the  total  surface  area  of  nuclei,  and 
accelerate  heterogeneous  nucleation,  resulting  in  a higher  nucleation  rate  at 
a lower  super saturation i Larger  values  of  the  factor  also  mean  higher 
concentrations  of  adsorbed  monomers  on  the  substrate  and  sui  Increase  in  the 
nucleation  rate. 

The  values  of  v(Rjj)  and  ^d(Pd)  somewhat  xincertain  euid  depend 
largely  vpon  the  nucleation  situation,  but  fortunately  they  are  not  significant 
qviantltles  eis  seen  in  these  figures.  The  features  of  heterogeneous  nucleation 
and  condensation  are  dominated  largely  by  the  contact  angle,  6. 
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3.  COWDMSATION  OF  WATER  VAPOUR  ON  HETEROGENEOUS  NUCLEI 


The  approeu±  to  the  problem  of  condensation  of  water  velour  in  a 
shock  tube  is  the  sane  as  In  Ref.  1.  A mixture  of  water  vspour  and  air,  which 
contadns  solid  particles  acting  as  heterogeneous  nuclei,  is  sviddenly  expanded 
by  a rarefaction  wave  in  a shock  tube  and  cooled  \mtil  the  vapour  becomes 
s\per saturated.  Heterogeneous  and  homogeneous  nucleation  and  condensation 
proce^sses  start.  The  associated  variations  of  pressure,  tenperature,  siper- 
saturatlon,  condensate  mass  and  other  physical  qiiantities  can  be  obtained  by 
solving  the  gasdynanlc  equations  covpled  with  the  nucleation  rate  equations 
obtedned  in  the  preceding  chapter. 

3.1  Nimerical  Model 

A mixture  of  water  vapour  and  air  or  an  inert  carrier  gas  is  initially 
at  rest  in  the  high-pressure  driver  section  of  a shock  t\ibe.  At  time  t = 0,  the 
diaphragm  separating  the  driver  section  from  the  low-pressvire  channel  section  of 
the  shock  tube  is  suddenly  r\ptured,  and  the  mixture  is  expanded  into  the  channel. 
When  sufficient  siper saturation  occurs,  embryos  of  water- droplets  nucleate  heter- 
ogeneously on  the  substrate  of  solid  particles  contedned  in  the  carrier  gas  and 
then  grow  to  larger  water  drops  by  condensation  of  water  vapour.  Concerning  the 
gasdynamlc  and  thermodynamic  features  of  this  system,  the  following  assimptlons 
are  also  made: 

(i)  The  effects  of  molecular  transport  leading  to  viscosity,  heat  conduction 
and  diffusion  are  neglected. 

(li)  Solid  particles  move  at  the  same  speed  as  the  mixture  regardless  of 
water-vapour  condensation  on  them. 

(ill)  The  thermodynamic  properties  of  the  mixture  are  the  weighted  sums  of  the 
corresponding  properties  for  the  single  system. 

3.1.1  Thermodynaaiic  Properties 

Let  be  the  density  of  species  a,  and  pa*  its  molecular  weight. 

Denote  the  subscripts  v,  i,  s and  i as  vapour,  liquid  phase  of  water,  solid 
particle  and  inert  carrier  gas,  respectively.  The  equation  of  state  for  a 
mixture  of  water  vapour  and  carrier  gas  is. 


where  ft  is  the  universal  gas  constant.  The  mixture  density  is  written  as, 

P = Pi  + Pv  + P]  + Pg  (^2) 

where,  and  Pg  are  the  masses  of  liquid  phase  of  water  and  solid  particle  per 
unit  vdvime  of  the  mixture,  respectively.  By  defining  the  initial  specific 
humidity  Uq,  the  condensate  mass  freu:tion  g and  the  mass  fraction  of  solid 
particle  gg  as. 


u 

o 


.%* '‘I 


0*3) 


the  equation  of  state,  Eq..  4l,  can  be  rewritten  as, 

g 

i +-il 

' '1  *“v 

The  specific  enthalpy  of  the  nixt\are  is. 


(i)  - g \ 

2+-^  )pftl 


'’l  '’v  H ‘’a 

+ TT  4,  + h,  + — h 

pipvpips 


(44) 


(45) 


where,  ha  is  the  specific  enthalpy  of  a-phase.  By  using  the  first  law  of  thermo- 
dynamics. 


\ - L + c^(T^  - T)  ~ h^  - L 


where,  L is  the  latent  heat  of  vapourization  and  by  taking  the  water  vapour  and 
the  carrier  gas  to  be  thermodynamically  perfect  gases,  the  specific  enthalpy 
and  the  i^eclflc  Internal  energy  of  the  mixture  can  be  expressed  as. 


h = Cpo  T - gL 


e = T - gL  - £ 


where 


°D  ® (l  - W - g )Co.  + W C + ff  ^ 

Bo  ' o ^i  o Pv  *s  T 

The  frozen  sound  speed  in  an  ideal  gas  mixture  is  defined  as, 


A . 1^81.'^  1 

l\ 

Using  Eq.  46,  the  sound  qpeed  is  then  e3q>ressed  as. 


(46) 


(47) 


(48) 
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where, 


, 1 - CJ  - g u 

1 o s ^ _o 


The  sv5>ersaturation  s defined  as  the  ratio  of  the  pressure  of  the  vapour 
phase  to  the  saturation  pressure  is  given  by, 


P W - 

s =3  _2__I 

where,  pg  for  water  vapovir  is  (Ref.  l). 


(^9) 


p = 10 

■^s 


(b-a/t) 


A = 2263 .0°K  B = 6.064 


(50) 


3.1.2  GasdynaMc  Equations 

By  neglecting  the  effects  of  molecular  transports  according  to  assuap- 
tion  (i),  the  equations  of  continuity,  momentum  and  energy  are  expressed  as 
follows : 


^•c|j(pu)  =0 


(51) 


I- 


du  . du 

3t  ^ "35? 


de  d ( 1 
dt  ^ dt  V p ^ ~ 


where,  u is  the  velocity  of  the  mixture  and 


dt 


1 ap 

P ^ 


3 3 

3t 


(52) 


(53) 


These  differential  equations  covgpled  with  the  nucleation  rate  equations 
can  he  solved  numerically  using  the  method  of  characteristics.  As  stated  in 
Ref.  1,  however,  this  method  is  not  valid  for  flows  when  the  characteristics 
intersect  to  form  shock  waves.  For  handling  such  a problem.  Lax  method  can  be 
used  by  smoothing  the  discontinuity  through  the  nature  of  its  differencing  scheme 
by  inplicitly  introducing  a viscosity  into  the  equations  of  motion. 

Lax’s  scheme  reqviires  all  differentisil  equations  to  be  expressed  in  the 
Lagrangian  form,  in  which  the  coordinates  and  state  variables  of  a flow  particle 
are  a function  of  time  t and  a parameter  6,  which  identifies  the  particle.  By 


J 
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taking  the  spatial  coordinate  at  the  tiiae  origin  t » o,  as  the  parameter,  the 
equation  of  continuity  can  be  written  as. 


(l»0)  - 

(TO 


ST 


x(l,t) 


(54) 


where,  p(|,t)  is  the  density  of  particle  | at  time  t and  x(|,t)  is  the  coordinate 
of  its  spatial  position  [I  = x(S,0)].  Since  the  particle  velocity  is, 


u=3^x(6,t) 


Equation  54  can  be  rewritten  as. 


(55) 


Using  Eqs,  54  to  56,  the  momentvim  (52)  and  energy  ^3)  equations  reduce  to  ’ 


3u 

ST 


1 

p(t;ot 


(57) 


Se  . p . ^ 

5t 


(58) 


The  last  equation  is  further  rewritten  by  using  Eq.  57  as, 


( e + I (pu)  (5841) 

Equations  56  to  58  together  with  the  rate  eq\aations,  Eqs.  38  to  40, 
and  the  thermodynamic  eqviations,  Eqs.  44  to  50,  describe  the  noneq\3ilibrium 
heterogeneous  condensation  of  water  vapour  in  one-dimensional  nonstationary 
flows.  They  can  be  solved  numerically  using  Lax  method  following  the  same 
technique  used  in  Ref.  1. 

3.1.3  Nianerical  Model 


Following  Ref.  1,  nondimensional  quantities  are  now  introduced, 


P 


f =1. 

"o’ 


u = 


u 


T = — 


y . p(y»°)  I, 

y Pc 


(59) 


PJ  A,  8 . p,  B, 


0 * P«  r' 


C,  <R^  = <R^ 
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where,  quantities  with  subscript  o refer  to  the  condition  in  the  driver  section, 
is  the  speed  of  soxmd  and 


where 


i = a t , 
o r’ 


r =a  1 p. 


The  gasdynamic  equations  become; 


St  5y 


(60) 


St  Sy  ^o 


(61) 


where 


M ^ 

St  by 


V = -,  E = e + i u^,  7=a^~2 

2 ’ 'o  o p 


1 
_ 3 

P 


The  thermodynamic  equations  are, 

H. 


E 


-iu-{ 


iT-®  js;* 


(r„  - 1)  ^ g - jr4r  e 

V -^'o'  o ' 'o 


1^0  - 

p = p(to  _ g)  __  T 
■^v  '■  o u 

*^v 


= TpAi(l-l/T) 


P 10 

o V 


Ai  = 2263 .O/T 


(62) 


s = 


(63) 

(64) 

(65) 

(66) 
(67) 


X = 


Po  o 


o Pso  ° 
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The  rate  equations  can  be  written  as, 


k = < 1^2  > ^ 


at 

= kir  Q ■> 
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• • 
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The  surface  tension  is  assumed  to  vary  linearly  with  ten5)erature  in  the 
same  way  as  used  in  Ref.  1, 


0-  = — = 1.693  - 0.00254T..  f 
r ® 


u = 75.6  dyne/cm^ 


The  latent  heat  of  condensation  is  taken  as  L = 623  cal/g.  Other  physical  quan- 
tities are  the  same  as  those  used  in  Ref.  1,  and  given  in  Tables  1 and  2 of  the 
paper.  Concerning  the  physical  properties  of  solid  particles,  those  of  S-O 
are  used,  although  they  have  negligible  effect  on  the  overall  properties  of^the 
mixture. 

Detailed  data  on  contact  angles  in  the  system  of  solid  particles  of 
Aitken  nuclei  wetted  by  water  are  not  available.  The  angles  0 = 60*,  90°  and 
120°  are  taken  as  representative  values.  Althovigh  the  values  of  the  vibrational 
frequency  of  adsorbed  atoms  and  the  free  energy  of  desorption  are  not  certain, 
they  are  roughly  estimated  as  v = (O.l  ~ 10)  10^®/ sec  and  = (2  ~ 200)  10® 
cal/g-atom  (Ref.  8),  that  is, 


Rjj  = 10“*  ~ 10®  cm,  ^ .r  3 ~ 300 


Equations  60  to  74,  with  the  foregoing  physical  properties,  were  solved 
manerically  for  the  same  case  used  in  Ref.  1 by  enploying  the  same  Lax  conputa- 
tlonal  technique.  However,  in  the  present  cslculations , the  mesh  sizes  of  spewre 
and  time  coordinates  were  changed.  The  size  of  space  mesh  chosen,  <iy  = l/4O0, 
which  corresponds  to  a distance  of  about  1 nm  between  two  particle  paths  (cf. 

0.5  nm  in  Ref.  l) , Doubling  the  size  of  space  mesh  is  Associated  with  the 


19 


trxmcation  error  of  about  5 ~ con^ared  with  the  results  of  I/8OO  mesh  size, 
but  it  decreases  the  con?)utational  time  by  a factor  of  about  three  or  four.  The 
stability  condition  time- increment  constant. 


I = ^ (p  a.) 
c ^ fW 


was  selected  as  - 0.25  (cf.  0.75  in  Ref,  l)  in  order  to  diminish  the  error 
associated  with  the  increased  space-mesh  size,  A Lax  dissipation  factor, 

M =0.2,  was  used.  The  numerical  conputations  were  performed  on  an  IBM  370-165 
at  the  Institute  of  Conpxiter  Science,  University  of  Toronto. 

3»2  Numerical  Results  and  Discussions 

3*2.1  Cccpaxison  with  Homogeneous  Nucleation 

Figures  11  to  15  show  the  resvilts  of  the  present  munerical  calcvilations 
of  heterogeneo\is  condensation  of  water  vapour  in  a shock  tube.  They  axe  conpaxed 
with  those  of  homogeneous  condensation.  The  numerical  constants  of  heterogeneous 
nucleation  axe 


< R^  > = 3.12  X 10”®  cm^  cm“®  (Curve  B in  Fig.  l) 

e = 90“,  Rjj  = < r2  >“1  = 0.311  X 10®  cm,  = 34 

Althoiigh  the  effects  of  these  constants  will  be  discussed  in  detail  in  the  next 
section,  the  obtained  results  axe  considerably  Influenced  by  the  value  of  the 
contact  angle. 

The  predicted  flow  field  of  the  nonequilibrium  nonstationaxy  rarefaction 
wave  in  the  physical  (x,  t)-plane  is  shown  in  Fig.  11.  The  head  of  the  rarefac- 
tion wave  moves  into  the  mixture  at  rest  with  the  state  sound  speed.  Behind  the 
rarefaction,  the  mixture  is  cooled  to  svpersatviration.  Consequently,  nucleation 
takes  place  on  the  heterogeneous  nuclei  with  an  appropriate  time  lag.  In  the 
figure,  the  condensation  wave  is  defined  as  the  locus  of  points  along  psurticle 
paths  vdiere  the  sipersaturation  reaches  its  maximum v alues  (cf.  Fig.  l4b). 

Owing  to  the  heat  release  associated  with  condensation,  the  condensation  wave 
is  followed  by  a shock  wave,  which  is  strictly  defined  by  the  locus  of  intersec- 
tion points  of  left-running  characteristics  (see  Ref.  l).  In  the  present  represen- 
tation, however,  the  shock  wave  is  defined  by  the  locus  of  points  of  the  pressure 
maximum  after  svper cooling.  Between  the  condensation  wave  front  and  the  shock 
wave,  there  exists  a condensation  zone,  hatched  in  the  figure,  where  the  flow 
variables  change  appreciably  due  to  heat  release.  It  can  be  seen  from  the  fig\are 
that  condensation  in  the  heterogeneous  case  (0  = 90°)  occurs  sooner  than  in  the 
homogeneous  case  (HM)  . The  width  between  two  fronts  is  also  narrower  for  0 = 90*. 
The  frozen  isentroplc  tall  of  the  wave  corresponds  to  a frozen  Mach  number  for  a 
diaphragm  pressure  ratio  (li^^  = 0.784), 

Figures  1^  to  15  show  time  and  space  variations  of  the  pressure,  tenw 
perature,  s\p)ersaturation  and  condensate  mass  fraction,  a)  along  three  particle 
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paths.  Initially  at  -5  (a),  -10  (b)  and  -20  (c)  cm  from  the  diaphragm  of  the 
shock  tube,  and  b)  at  four  time  lefvels,  0.48  (a),  0.72  (b),  O.96  (c)  and  1.2 
(d)  msec.  As  shown  in  Figs.  12  to  13,  in  heterogeneous  condensation,  the 
decreaise  in  the  pressure  and  tenperature  before  the  onset  of  condensation  is 
smaller  than  that  in  homogeneous  condensation.  The  successive  rise  of  their 
values  after  the  onset  of  condensation  is  also  smaller . The  onset  of  condensa- 
tion occurs  for  particle  paths  farther  from  the  diaphragm  or  at  lower  time, 
levels.  Space  and  time  variations  of  super saturation  are  shown  in  Fig.  14.  It 
is  readily  seen  that  heterogeneovis  nucleation  leads  to  lower  values  of  maximum 
sxqpersaturatlon  since  the  onset  of  condensation  occurs  closer  to  the  wave'head. 
Variations  of  the  condensate  mass  fraction  are  plotted  in  Fig.  15.  In  the 
heterogeneoiis  cases,  the  condensate  mass  fraction  increases  more  sharply  15)  to 
the  equilibrium  value. 

Since  the  svqjer saturation  is  directly  related  to  the  nucleation  rate, 
the  latter  attains  its  maximum  approximately  at  the  highest  siqier saturation. 
Thereafter,  the  condensation  of  water  vapour  onto  the  solid  nuclei  plays  an 
increasingly  inportant  role  in  the  phase- changing  process.  The  increased  heat 
release  due  to  condensation  causes  an  increase  in  the  tenperatvjre  of  the  mixture. 

In  heterogeneous  nucleation,  the  activation  energy  is  greatly  reduced  by  the 
fact  that  the  interfacial  free  energy  of  the  chemical  bond  between  vapovtr  atoms 
and  substrate  is  larger  than  that  between  liquid  atoms  and  substrate  (t^gy  > a^) 
and  their  difference  reduces  the  overall  interfacial  energy  for  formation  of  the 
liquid- vapoiu-  interface  of  enibryos.  This  acts  equivalently  to  decrease  the  surface 
tension.  Less  activation  energy  restilts  in  a higher  nucleation  rate  with  lower 
sirpersaturation,  lower  stper cooling,  a higher  propagation  velocity  of  the  conden- 
sation wave,  and  a narrower  condensation  zone. 

3«2.2  Effects  of  Physical  Properties  of  Nuclei  ' 

It  is  shown  in  Chapter  2 that  the  size  distrib\ition  of  condensation  j 

nuclei  hardly  affects  the  condensation  process  in  the  expansion  of  the  mixture,  ! 
U3iless  their  size  is  coinparable  to  that  of  nucleating  embryos.  For  such  a j 

size  distribution  as  that  of  Aitken  nuclei  shown  in  Fig.  1,  the  property  which  j 

contributes  directly  to  the  nucleation  process  is  the  total  surface  area  of  the  | 

nucleus  < R®  >.  Other  factors  controlling  the  heterogeneovis  nucleation  are  the  ] 

concentration  of  adsorbed  monomers  on  the  surface  of  the  substrate  and  the  contact  I 
angle  of  embryo.  The  former  is  characterized  by  the  vibrational  frequency  of  ] 

adatoms  v(Rh)  and  the  desoiption  free  energy  of  adsorbed  monomers  ^d(Pd)*  Since 
the  activation  energy  of  nucleation  is  much  larger  than  the  desorption  energy  of  i 
adsorbed  monomers,  the  latter  does  not  greatly  affect  the  nucleation  process  but 
reduces  slightly  the  overall  activation  energy  of  nucleation.  Thus,  the  factor  j 
RH  has  a physical  meaning  related  to  the  concentration  of  monomers  adsorbed  onto  i 
the  svirface  of  the  substrate.  It  is  a constant  coefficient  of  the  nucleation 
rate  n.  For  the  sake  of  reducing  conputer  costs,  only  the  contact  angle  0,  and 
the  factor  Rjj,  were  considered  as  the  physical  parameters  that  affect  the  flow  1 
field  of  the  nonequilibrium  nonstationary  rarefaction  wave  of  water-vapour/carrier-j 
gas  mixtiore  in  a shock  tube.  In  the  expression  of  the  condensate  mass,  however,  j 
the  factor  Rji  is  always  covipled  with  the  overall  surface  area  < R®  >.  The  influena 
of  the  overall  surface  area  or  the  overall  ninnber  of  nuclei  can  be  inferred  by  | 
inspection  of  the  effect  of  the  factor  Rjj.  , 

Figures  I6  to  21  show  the  effects  of  the  contact  angle  0 and  the  factor  : 
Rh»  For  cases  to  conpare  the  effects  of  the  contact  angle,  Rg  = 0.3II  x 10^  cm 


is  used.  For  the  curves  denoted  by  Rg,  Rjj  ■ 0.311  x 10“*  an  and  9 = 90*  are 
used.  The  flow  fields  in  the  (x,  t)  plane  are  shown  in  Figs.  16  and  17,  the 
effects  of  the  contact  angle  (0  = 90°,  120®)  and  the  factor  Rh  (Rr  10®,  10** 
cm),  respectively.  For  larger  contact  angles,  which  means  embryos  of  n»re 
^hertcal  shape,  the  heterogeneous  nucleation  loses  its  influence  and  tends  to  have 
the  features  of  homogeneous  nucleation.  Reducing  the  value  of  the  factor  R„  leads 
to  a lower  nucleation  rate  and  higher  supercooling.  Figures  I8  to  21  are  the 
time  and  space  variations  of  the  pressure,  tenperature,  supersaturation  and  con- 
densate mass  fraction,  respectively.  In  the  curves  denoted  by  R«,  the  factor  R„ 

Is  of  two  Orders  less  tiian  'tha't  of  other  ctirves#  It  does^  however,  have  an 
influence  less  than  that  of  0 = 120® . It  is  the  contact  angle  that  really 
^aracterizes  the  significant  featvire  of  heterogeneous  condensation.  It  should 
be  noted  that,  for  an  accurate  prediction  of  the  flow  field  of  heterogeneous 
condensation,  the  contact  angle  of  enbryos  nucleating  on  the  substrate  of  solid 
particles  must  be  carefuUy  estimated.  Yet,  other  parameters,  including  the 
Derail  sitrface  area  of  the  nuclei  can  be  roughly  estimated  without  serious 
ejects.  Nevertheless,-  there  is  no  precise  knowledge  to  make  reasonably  good 
choices  possible  \mtil  existing  ejqperimental  data  from  various  experiments  are 
correlated. 

3.2.3  Conparison  with  E3<perim<»ntal  Data 

The  numerical  results  can  be  conpared  with  Kalra*  s experimental  study 
% ^ nOTstationary  expansion  of  water-vapour-nitrogen  mixture  in  a shock  tube 

"the  sane  considerations  of  the  pressure  history  discussed 
jf*  V*  space  location  of  the  experimental  measurement,  x = -17.0  cm  from 

the  diaphragm,  is  selected  as  corresponding  io  x = -23.4  cm.  To  obtain  a 
reasonable  fit  with  the  experimental  data  in  Ref.  1,  the  surface  tension  was 
+ ^alue  of  the  semi-enpirical  expression,  Eq.  75,  cr  = 82  dynes/cm 

to  00  dynes/cm,  with  the  conclusion  that  choosing  a suitable  value  for  the  surface 
tension  provides  a good  fit  for  the  experiments.  As  shown  in  Fig.  22,  in  hetero- 
geneous condensation,  a similar  fit  can  be  obtained  by  replacing  the  surface 

contact  angle.  By  choosing  a suitable  value  for  the  contact  angle 
^ experimental  data  can  be  fitted  fairly  well  except  for  the  stper- 

saturation  at  the  onset  of  condensation,  which  is  shown  in  brackets  in  Fig.  22, 

Figure  23  is  a conparison  with  the  experimental  data  as  plotted  for  the 
(^f  supercooling  and  the  cooling  rate  derived  by  Kalra 
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where,  Tg  is  the  tenperature  of  eqxiilibrium  saturation,  T.  the  tenperature  at  the 
onset  of  condensati^  and  (-dl/dt)  the  rate  of  cooling.  ^Concerning  this  relation. 

numerical  results  are  quite  appreciable.  The  best- fit 
conditions  of  the  pressure  history  in  Fig.  22,  0 = 90®,  in  the  heterogeneous  case, 

homogeneous  case,  result  in  the  poorest  agreement  in 
Fig.  23.  ^ng  to  the  limited  experimental  data,  it  is  not  certain  whether  this 
is  due  to  the  inadequacies  of  the  models  or  features  of  the  experimental  data, 
such  as  the  estimation  of  tenperature  from  the  pressure  profiles,  for  exaaple, 
or  whether  the  waves  are  planar.  It  cannot  be  concluded  at  present  which  process, 
heterogeneous  or  homogeneous  condensation,  adequately  explains  the  experimental 
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It  should  also  be  noted  that  the  best  fits  for  homogeneous  and 
heterogeneous  nucleation  and  condensation  will  give  points  very  close  to  each 
other  on  the  time-distance  plot  shown  in  Sisliem's  Fig.  ^9  (Ref.  1).  Conse- 
quently, the  schlieren  records  of  Glass  and  Patterson  (Ref.  U)  cannot  be  used 
as  a decisive  test  for  the  two  coopeting  theories,  since  the  traijectory  of 
the  condensation  shock  agrees  with  either  point  (only)  from  the  two  theories. 

An  inspection  of  Fig.  22  shows  that  the  experimental  time-vcu'iation 
of  pressure  at  the  onset  of  condensation  is  extremely  sharp  cooqpared  with  the 
numerical  results.  This  is  not  surprising  in  view  of  the  fact  that  the  inplicit 
artificial  viscosity  used  in  the  numerical  analysis  tends  to  spread  the  conden- 
sation and  shock  fronts.  It  is  also  worth  noting  that  factors  such  as  coagulation, 
solubility,  or  electrical  charges  have  not  been  treated  in  this  analysis.  Physical 
data  in  these  areas  are  lacking. 

4.  CONCLUDING  REMARKS 

By  using  a macroscopic  model  of  heterogeneous  nucleation,  a theoretical 
study  was  made  of  the  condensation  of  water-v^our/carrier-gas  mixtures  in  a 
nonequilibrivim  nonstationary  rarefaction  wave  generated  in  a shock  tube.  The 
present  analysis  is  conpared  with  the  results  of  homogeneous  nucleation.  Nuclea- 
tion is  assumed  to  take  place  heterogeneously  on  idealized  smooth,  spherical 
solid  particles  of  Aitken  nuclei,  which  are  chemically  and  electrically  inert. 

In  the  process  of  heterogeneous  nucleation,  the  controlling  factors 
are  the  size  distribution  of  nuclei,  the  concentration  of  monomers  adsorbed  on 
the  surface  of  the  substrate  and  the  contact  angle  of  embryos.  Of  these  factors, 
the  dominant  one  is  the  contact  angle.  By  decreasing  the  contact  aujgle,  the 
activation  energy  of  nucleation  is  greatly  reduced  due  to  the  fact  that  the  forma- 
tion of  a liquid-vapour  interface  is  replaced  by  a liquid- solid  interface  having 
a lower  interfacial  free  energy.  At  less  supercooled  states  and  lower  siper- 
saturations,  the  onset  of  condensation  occurs  closer  to  the  wave  head  where 
tenperatvtres  and  pressures  are  higher.  The  reverse  is  true  for  homogeneous 
condensation . 

When  nuclei  are  in  the  size  range  of  one  or  two  orders  larger  than 
that  of  the  embryos  (R  ~ 10“  ‘cm)  and  the  value  of  the  contact  angle  is  not  too 
large  (0  < 120“ ),  heterogeneous  nucleation  dominates  the  condensation  process 
and  homogeneous  nucleation  can  be  neglected.  In  this  case,  the  size  distribution 
of  nuclei  affects  the  overall  nucleation  and  condensation  only  through  the  total 
surface  area  of  the  nuclei.  The  concentration  of  monomers  adsorbed  on  the  surface 
of  the  substrate  is  dcminated  by  the  vibrational  frequency  and  the  desorption  energy 
of  adatoms.  These  act  as  a factor  on  •vdilch  the  nucleation  rate  has  a linear 
dependence.  The  reduction  in  the  total  siirface  area  of  nuclei  and  the  factor 
depending  on  the  surface  concentration  of  monomers  results  in  weakening  the  effect 
of  the  contact  angle. 

The  nucleation  rate  is  very  sensitive  to  the  contact  angle.  Consequently, 
its  value  should  be  carefully  chosen.  An  appropriate  value  for  the  contact  angle 
makes  it  possible  to  fit  the  numericeil  results  with  experimental  data.  From  the 
data  available  at  ITTIAS,  a contact  angle  of  90*  appears  reasonable.  However, 
this  does  not  substitute  for  an  actual  measurement.  Perhaps  it  is  only  indicative 
of  an  averaged  result,  0 = 90‘’»  This  is  such  an  inportant  quantity  for  hetero- 
geneous nucleation  that  it  makes  an  experimental  determination  inperatlve. 


In  the  present  theory,  the  shape  of  enibryo  and  substrate  different 
from  a sphere,  the  electric  charge  and  the  solxibility  of  nuclei  were  not  taken 
into  account.  These  are  inportant  factors  that  can  significantly  redvice  the 
activation  energy,  thereby  increasing  the  heterogeneous  nucleations.  To 
examine  the  effects  of  svich  factors,  the  present  theory  can  be  readily 
extended.  However,  other  uncertainties  wo\ild  be  introduced  when  selecting 
suitable  values  for  the  corresponding  physical  properties. 

When  nuclei  or  embryos  have  a higpier  number  density,  coagulation 
owing  to  Brownian  notion  takes  place  and  is  another  factor  to  be  considered. 

By  assuming  the  coagulation  constant  as  10“^  cm®/(sec  particles)  (Ref.  12)  and 
the  characteristic  time  as  lO'®  sec,  the  number  density  above  which  coagulation 
has  a considerable  effect  can  be  estimated  as  10^®  parti cles/cm®.  Since  the 
assxaned  nvunber  density  of  heterogeneous  nuclei  is  about  10®  parti cles/cm®, 
coagulation  hardly  affects  the  results  within  such  a characteristic  time.  In 
the  homogeneous  nucleatlon  case,  however,  it  may  be  a factor  to  be  taken  into 
account  because  of  the  higher  number  density  of  nucleated  droplets  (>  10^ ^ 
parti cles/cm®,  see  Ref.  l).  Kalra's  experimental  results  can  be  explained 
equally  by  both  homogeneous  and  heterogeneous  nucleation  theories  with  suitable 
values  of  surface  tension  and  contact  angle,  respectively.  In  view  that  the 
former  requires  a lower  value  of  surface  tension  and  a possible  consideration 
of  coagulation,  heterogeneous  nucleation  and  condensation  wo\ild  be  closer  to 
the  experiments.  However,  the  fact  that  the  contact  angles  have  not  been 
measured  throws  some  doubt  on  that  analysis.  It  must  therefore  be  concluded 
that  the  choice  of  model,  heterogeneous  or  homogeneous  condensation,  can  only 
be  settled  by  a decisive  set  of  experiments. 
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TABLE  2 

INITIAL  CONDITIONS  IN  THE  SHOCK  TUBE 


Driver  Conditions 

MLxtiire  of  water  vapour  and  nitrogen 

Wjj  = 0.0176 

<Pq  = 0.9727 

= 18.6  mn  Hg 
Pgo  = 19.23  imn  Hg 
= 680  mn  Hg 

= 0.1025  X 10“®  g/can® 

T^j  = 295 .3"K 

a^  = 352.2  m/sec 

Channel  Conditions 

Air 

px  = 100  mn  Hg 
Ti  = 295 .3'K 


NucIm*  Radius  R cm 

FIG.  1 SIZE  DISTRIBUTION  OF  NUCLEUS  PARTICLES 


2 CAP-SHAPED  EMBRYO  NUCLEATING  ON  SURFACE  OF  INSOLUBLE  SUBSTRATE 
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FIG.  4 COEFFICIENTS  (f.  K)  IN  NUCLEATION  RATE  EQUATIONS 
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FIG.  15-a  TIME-VARIATION  OF  CONDENSATE  MASS  FRACTION  ALONG 
PARTICLE  PATHS  INITIALLY  AT  -5(a),  -10(b),  -20(c).  cm  THE  DIAPHRAGM 


FIG.  15-b  SPACE-VARIATION  OF  CONDENSATE  MASS  FRACTION  FOR 
TIME  LEVaS:  0.48(a).  0.72(b).  0.96(c).  1.2(d).  msec. 
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FIG.  16  CONDENSATION-WAVE  STRUCTURE;  EFFECT  OF  THE  CONTACT  ANGLE  e 

(Rj^  - 0.311  X 10®  cm) 


\Mbvt  Front 


y^lsantropic 
\ Rarvtaetk 


Rarvtaetion 
Wovt  Tail 


CondansQtion. 
Wove  Front 


1h  ■0.311x10^  cm 


RH*0i3llxlO*em^ 


Rorafaction  \Ma««  Hiad< 


FIG.  17  CONDENSATION-WAVE  STRACTURE;  EFFECT  OF  THE  FACTOR  R^^  (e-  90*) 


FIG.  18-a  TIME-VARIATION  OF  PRESSURE  ALOTJG  PARTICLE  PATHS 
INITIALLY  AT  -5(a),  -10(b),  -20(c)  cm  FROM  THE  DIAPHRAGM 
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FI6.  22  COMPARISON  WITH  EXPERIMENTAL  DATA; 
TIME-VARIATION  OF  PRESSURE  AT  STATION  -23.4  cm  FROM  THE  DIAPHRAGM 
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FIG.  23  COMPARISON  WITH  EXPERIMENTAL  DATA;  SUPERCOOLING 
VERSUS  COOLING  RATE 


Referring  to  Fig.  2,  one  obtains  the  geometriceJ.  properties  of 
shaped  embryo  on  the  spherical  surface  of  a solid  particle  as  follows, 


Liqiaid-vapour  surface  area: 


= hv  b 

V V 

Liquid- solid  surface  area: 


where 


The  coefficients  a,  by  and  bg  are  shown  in  Fig.  3 as  a function  of  r/R. 
r/R  < 0.1,  these  coefficients  take  on  almost  constant  values, 

a^  = (l  - COS0)®  ^ 7 ^ 

\co  = I 

^Soo  “ ¥ 


a cap‘ 


(A- 


(A- 


(A- 


For 


(A- 


APPENDIX  B:  INTEGRATION  OF  RATE  EQUATIONS 


The  rate  of  production  of  the  liquid  phase  by  heterogeneous  nucleation 
and  condensation  is  given  by  Eqs.  3^,  • 


f 

dt  J 


3T 


(B-1) 


sr  ■ { f { I '*0  - It  ( 

o 

where,  the  svirface  area  of  a nucleus,  wetted  by  condensate  droplets,  Ag  is  given 

By, 

^ ■"/{  "^^o  - 3T  \ 

o 

The  liquid-solid  interfacial  area  of  a condensate  droplet  nucleated  at  time  t is 
given  by, 

a^(T,t)  = a*(r)  + ag(r)dT' 

T 

where,  ag  is  the  growth  rate  of  the  liquid- solid  surface  area. 

For  the  numerical  integration  of  these  eq\iations,  it  is  possible  to 
reduce  them  to  a set  of  differential  equations.  Successive  differentiation  of 
Eq.  B-2  with  respect  to  time  leads  to, 

3F  = '"■  fi!  { (A„  - A,)  f a . A ]• 

i ■ I (A„  - A^r*=  b,  . B 

^=-(A-A)r*c.+rC 
dt  p ' o s'  D o 

(B-5) 

i = ? - Ag)<=c  "o 

= f (^o  - ^ 


'tJ.  . •TJ^'«  4?<B7C51> 


HlCVJ 


where, 


o 

o 

C s I 2 (A„  - A„)  c„  dT 
J p o s c 

o 

Cl  s / S (A^  - A ) c dT 
J p ' O s'  Cl 

o 

•/  f <=02 


(B-6) 


where 


r c®  ', 

= K -j^  1 - Cg(cos0  - j)  * ^ ~ C08®e)j- 


S(rc- ) 

C = K -A 
c or 

^c 

c = K 
Cl  dr 

dc 

%2  sr 


(B-7) 


"7  V ^3S  3?; 

o 

{i  - Oj  (cose  - |)j-[  ji  . 0,  (cose  - 1)|  +^(ccse  -|)j. 


0 1-1 
+ ^ ^ COS0(l  - COS0)(l  - cos^0) 


(B-8) 


B-2 


J 


Vhen  hardly  changes  with  respect  to  r or  r/R, 

=nH.l  -J*f  (*o  - ^ = 0 


which  leads  to 


-Sr  ” f ^^o  - °cn 

The  integral  equation  can  be  replaced  then  by  sequential  ordinary  differential 
equations  in  a closed  form. 

Similarly,  the  Integral  equation  of  the  wetted  surface  area  of  a 
nucleus  (B-3)  can  be  reduced  to  the  following  sequence  of  ordinary  differential 

equations. 


'rr^  = n(A  - A ) kirr*^  b + r B 

Ot  ' O s'  SOS 

3b  . 

- A ) 2ttt*  c . + r C 
ot  ' o s'  sb  o s 


35“  * n(A  - A ) 2v  o + r C 
ot  ' o s'  sc  o Si 


-Ttft  =•  n(A_  - A ) 2tt  c + r C 
ot  ' o a SCI  o s2 


= n(A^  ~ A ) 2v  c +r  C 
ot  ' o s'  see  o s 


“ I - a ) r c . dT 
s J ' o s'  sb 

C s I n(A„  - A ) c dT 
s J o s'  sc 

o 

C « I n(A_  - A ) c dT 
sx  J ' o s'  sex 

o 

pt 


C s 

S2 


C dT 
SC2 


f 


When  c takes  on  almost  constant  values  with  respect  to  r or  r/R, 
sen  ' 


which  yields > 


! -3T  = °(^o  ■ ^sn 

\ 

\ 

' It  means  that  the  integral  equation  (B-3)  can  be  explicitly  replaced  by  a set  of 

' ordinary  differential  equations. 


In  Fig.  5 are  shown  the  coefficients  K,  c^,  cc,  Cg-jj  and  Cg^  which  are 
approximately  constant  for  r/R  <0.1,  taking  the  following  values; 


In  this  case,  the  differential  equations  can  be  written  as  follows: 


% 

dt 

r*3 

3 

dA 

dt 

• 

_ n 
~ P 
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• A )r*2  b 
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+ r B 
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dB 
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_ n 
' P 
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r c 
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dC 
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" P 

(1  - 

■ 'o- 

dA 

= n(l  - A ) 47rr*®  b + ^ 5 
d.t  ' s soo  o s 

dB  . 

■37-  = n(l  - A ) 27rr*  c , + f C 

dt  ' s'  sboo  o s 

dC 

~ = n(l  - A^)  27T  c^^ 

dt  S SCoo 


where,  gp.  A,  B, ...  mean  g^.  A,  B,...  divided  by  Aq  which  is  equel  to  kHR^.  The 
rate  of  production  of  overall  condensate  mass  fraction  is  then  given  by. 
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